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Abstract — This paper develops a sum-power minimized re- 
source allocation (RA) algorithm subject to a sum-rate con- 
straint for cooperative orthogonal frequency division modulation 
(OFDM) transmission with subcarrier-pair based opportunistic 
decode-and-forward (DF) relaying. The improved DF protocol 
first proposed in [lj is used with optimized subcarrier pairing. 
Instrumental to the RA algorithm design is appropriate definition 
of variables to represent source/relay power allocation, subcarrier 
pairing and transmission-mode selection elegantly, so that after 
continuous relaxation, the dual method and the Hungarian algo- 
rithm can be used to find an (at least approximately) optimum RA 
with polynomial complexity. Moreover, the bisection method is 
used to speed up the search of the optimum Lagrange multiplier 
for the dual method. Numerical results are shown to illustrate 
the power-reduction benefit of the improved DF protocol with 
optimized subcarrier pairing. 

Index Terms — Cooperative communication, resource alloca- 
tion, decode and forward, OFDM. 



I. Introduction 

Cooperative orthogonal frequency division modulation 
(OFDM) transmission with subcarrier-pair based decode-and- 
forward (DF) relaying and associated resource allocation (RA) 
were studied in |2)-|[5] when the source-to-destination link 
exists. In Q-JH, an "always-relaying" DF protocol was 
considered. To better exploit the frequency-selective channels, 
we have proposed an opportunistic DF relaying protocol in (Q, 
i.e, a subcarrier in the first time slot can either be paired with 
a subcarrier in the second slot for the relay-aided transmis- 
sion, or used for the direct source-to-destination transmission 
without the relay's assistance. A major drawback for that DF 
protocol is that, a subcarrier unused for relaying in the second 
slot becomes idle, which wastes spectrum resource. To address 
this issue, we first proposed in |fl~) an improved DF protocol, 
which allows the source to make the direct transmission over 
the subcarriers unused for relaying in the second slot. This 
protocol and its RA were later intensively investigated, e.g., 
in ll6l- lfTTl . Note that in (TJ, O — (9J, a priori subcarrier pairing 
was used, i.e., when the relay-aided transmission is used for 
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a subcarrier in the first slot, the same subcarrier in the second 
slot is always paired with this subcarrier for DF relaying. In 
(0-11), iflOl , ifTTl . the optimization of subcarrier pairing was 
considered. 

Recently, energy-efficient communication is becoming in- 
creasingly important |[T2l . In view of the fact that many 
existing works focus on spectral-efficiency maximized RA, 
we develop an RA algorithm for minimizing the sum power 
subject to a sum-rate constraint for the improved DF protocol 
with optimized subcarrier pairing (OSP). Compared with the 
algorithms designed in J4], ifTUl . our algorithm uses a new 
method to define indicator variables for representing subcarrier 
pairing and transmission-mode selection, by regarding the 
subcarriers for the direct transmission in the two time slots 
as virtual subcarrier pairs. Moreover, the bisection method is 
used to find the optimum Lagrange multiplier, which is faster 
than the incremental-update based subgradient method used in 

a, ed. 

Notations: C(x) = ^log 2 (l +x) and [x] + = max{x, 0}. 

II. System and protocol description 

Consider the scenario where a relay assists a source's 
transmission to a destination. The improved DF protocol in 
JT] is used. Specifically, every data-transmission session takes 
two consecutive equal-duration time slots and OFDM with K 
subcarriers is used. To facilitate description, a subcarrier used 
in the first slot is denoted by subcarrier k and one in the second 
slot by subcarrier I hereafter. In the first time slot, the source 
radiates OFDM symbols, using P s .k.i as the transmit power 
for subcarrier k. The source-to-relay and source-to-destination 
baseband-channel coefficients for subcarrier k are h ST) k and 
hsd,k, respectively. In the second slot, both the source and 
the relay synchronously radiate OFDM symbols, using P s ,z,2 
and P T ,i,2 as the transmit powers for subcarrier I, respectively. 
The relay-to-destination baseband-channel coefficient is h r d,i 
for subcarrier I. 

A subcarrier in the first slot can either be paired with 
one in the second slot for the relay-aided transmission, or 
be used for the direct transmission without relaying. Every 
unpaired subcarrier in the second slot is used for the direct 
transmission. In particular, if subcarrier I is used for the 
direct transmission, P Si ;.2 > is used while P r .z,2 = is 
imposed. The maximum average data rates over subcarriers k 
and I used for the direct transmission are C(P s .k.iG s d.k) and 
C(P s ,i,2G s d,i) bits/OFDM-symbol (bpos), respectively, where 



G 



and 



the noise variance for each 



subcarrier at every node's receiver. When subcarrier k is paired 
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with subcarrier I for the relay-aided transmission, the DF 
relaying is used in which case P s ,i,2 = is imposed while 
P&,k,i > and P r ,i,2 > are used (more details are available 
in [I]). Suppose P s ,fe,i + Pi,i.2 = P, it can readily be shown 
that the maximum data rate is equal to R k j = C(G k iP) bpos, 
where 



of variables (COV)_from P to P = {P k i, a kl , #w|Vfc, I}, 
where P ku a M and /3 M satisfy P u = ijyPfcj, a k i = t kl a k[ and 
f3 k i = tfciPki, respectively, V k, I. After collecting all variables 
into X = {I, P}, the RA problem can be rewritten as the 
problem (P2): 



if min{G aIik ,G ldt i} > G. 



d.k j 



G S r,fc 

is achieved when 



rmn{G sl% fc,G s d,A;} if min{G sr , fe , G ld j} < G sd ,k, 

(TJ. This maximum rate 



l^ar.fcl 



and G rd « = 



fc,i 



p 



Assume there exists a central control unit which knows pre- 
cisely {Gsr.fc, G s d,fc|V k} and {G r d.;|V I}, and determines the 
optimum RA (i.e., the source/relay power allocation, subcarrier 
pairing and transmission mode selection ) to minimize the sum 
power subject to the constraint that the sum data rate is not 
smaller than prescribed R lcq bpos. 

III. RA ALGORITHM DESIGN 

For any subcarrier assignment used by the improved DF 
protocol, suppose m subcarrier pairs are assigned to the relay- 
aided transmission, then it is always possible to one-to-one 
associate the unpaired subcarriers in the two slots to form 
K — m virtual subcarrier pairs for the direct transmission. 
Motivated by this observation, the RA problem is formulated 
by defining: 

. t\\ G {0,1} and P kl > 0, V k,l. tf t = 1 indicates that 
subcarrier k is paired with subcarrier I for the relay-aided 
transmission. When = 1, P k i is used as the total power 
for the subcarrier pair (k, I). 

. ig S {0,1}, a k i > and fi u > 0, V k,l. ijg = 1 
indicates that subcarriers k and I form a virtual subcarrier 
pair for the direct transmission. When = 1, P s . k .i and 
P s ,i,2 take the value of a k i and (3 k i, respectively. 

Let us collect all indicator and power variables in the sets I 
and P, respectively, and define S = {I,P}. The RA problem 
can be formulated as the problem (PI): 

m s in E(^ Pm + f ^ aki + tufa), 

k,l 

s.t. e {o,i},v k,l] 

E(*S + *S)=l.Vfc;X;(*«+*S)=l.Vi; 

I k 
Pki >0,a kl >Q 7 f3 kl >0,Vfc,i; 
/(S) > R Icqi 
where /(S) represents the maximum sum rate as 

/(S) = J2 {tfiC(G kl P kl ) + t%C(G sd , k a kl ) + t»C(G sd j(3 kl )) 

k,l 

Obviously, (PI) is a nonconvex mixed-integer nonlinear pro- 
gram. To find the optimum S, we first relax all indicator vari- 
ables to be continuous within [0, 1]. Then, we make the change 



mm 
x 



in P(X) = J2(Pki+aki+M 



k.l 



P if min{G srife ,G rd ,/} > G sd ,fe, 
if min{G sl! fc, G r d,/} < G sd>k , 



s.t. t»,ig e [0,l],V*,i; 

E(*S + *S) = i-va; E(*S + *S) = 1 ' VZ ; 

; k 
Pki>0,a kl >0,Ai >0,Vfc,/; 

- g(X) < --Rreq, 

where <?(X) represents the maximum sum rate expressed as 
g(X) = J2{4>(^i,Pki,G kl ) 

k,l 

+ 4>{^u ®ki, G s d,fe) + <i>(t kl , p k i, G sd) i)), 



and 



(1) 



Obviously (P2) is a relaxation of (PI). We will find an 
(at least approximately) optimum solution for (P2), and show 
that the S corresponding to this solution is still feasible, and 
hence (at least approximately) optimum for (PI). To this end, 
note that <f>(t, x, G) with fixed G is a continuous and concave 
function of t > and x, because it is a perspective function 
of C(Gx) which is concave of x fPH . As a result, <?(X) is a 
concave function of X in its feasible domain for (P2). This 
means that (P2) is a convex optimization problem. As can 
be checked, it also satisfies the Slater constraint qualification, 
therefore it has zero duality gap, which justifies the applicabil- 
ity of the dual method to find the globally optimum for (P2), 
denoted as X* hereafter. 

To use the dual method, /j, is introduced as a Lagrange 
multiplier for the rate constraint. The Lagrange relaxation 
problem for (P2) is the problem (P3): 



min L(n, X) = P(X) + /i ( P roq - g(X.) 
s.t. g [0,l],Vfc,Z; 



£(ts+*a)=i,v*;5[;(*£+*5)=i.v*i 

Pki >0,a kl >0,Ai >0,Vfc,/; 

where L(fi, X) is the Lagrangian of (P2). A global optimum 
for (P3) is denoted by X^ and the dual function is defined as 
d(fi) = L(/x,X /J ). Note that d(fi) is concave of fi > 0, and 
Prcq — 9(X- fl ) is a subgradient of i.e., V \i! ', d(fi') < 

d{(j) + (fj,' — fi)(R rcq — .g(X M )). The dual problem is to find 
the dual optimum /i* = argmin^o d(fx). 

Since (P2) has zero duality gap, two important properties 
should be noted. One is that fj,* > 0. This is because fi* 
represents the sensitivity of the optimum objective value for 
(P2) with respect to P rcq , i.e., — '- = fx H 1 3 1 - Obviously, 
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P(X*) is strictly increasing of R lcq , meaning that /a* > 0. 1 03). After knowing {££j|V k,l}, the optimum I can be 

The other is that fjb = /Lt* and X M = X*, if and only if X M is constructed according to the^way mentioned earlier. Finally, 
feasible and /Lt(<?(X M ) — R rcq ) = according to Proposition 
5.1.5 in |[T4l . Based on the above property, the \i > and 



X M that satisfies ,g(X p 



i? loq can be found as /i* and 



X*. Therefore, the key to using the dual method consists of 
two procedures to finding X M and fx*, respectively. We first 
introduce the one to finding X M as follows. 

1) To find X M when fi > 0: the following strategy is used. 
First, the optimum P for (P3) with fixed I is found and denoted 
by Pi. Define Xi = {I, Pi}. Then we find the optimum I to 
maximize L(/i, Xi) subject to the constraints on I in (P3). Xi 
corresponding to this optimum I can be taken as X M . 

Suppose I is fixed, it can readily be shown that the optimum 
P k i, a H and f3 H for (P3) are 

Ki = t^A((i, G kl );a kl = tgA(/i, G sd , fc ); Pki = *«A(m, G scM ) 



where A(/z, G) 



' k-i 

lo &2 e ,, _ i 
2 r G 



Using these formulas, 



Xi = {I, Pi} can be found. It can readily be shown that 

L(ji, Xi) = M i? req + {tfi A ki + t%B kl ) (2) 



k.l 



where 



A kl =A(/i, G kl ) - ii • C{G kl k(», G kl )) 
B kl =A{fi, G s d.fc) - n ■ C(G sd , fc A(>, G sd)fc )) 
+ A(/i, G scM ) - n ■ C(G scM A(/i, G scM ))- 

Now, it can be readily shown that the optimum I for (P3) 
is the solution to the problem (P4), 



I,{i/si|V k,l} 



,R ,D 



s.t. 6 [0,l],Vfc,Z; 

tki =ifci + ifcz, V fc, Z; 
y^Jkl = 1,V fc; y^tfc; = 



where extra variables {ifcz|V fc, /} are introduced. Note that 
ifcjA w + t^B M > t k iCki holds where C H — minjAfc;, B k {\. 
Let us label A k i as the metric for iS and as the metric for 
fj? z . This inequality is tightened when the entry in {^,ij? ; } 
with the smaller metric is assigned to t k i, while the other entry 
assigned to 0. This means that after the problem (P5): 

min y^tkiCki 
{tjtiVfe.J} *-f 

s.t. t H £ [0,l],Vfc,Z; 

y^Jm = 1,V fc; y^tfc; = 1, V f; 



the corresponding Xi = {I, Pi} is assigned to X p . Note that 
the Hungarian algorithm to solve (P5) has a complexity of 

0{K 3 ) na. 

2) To find fi*: an incremental-update based subgradient 
method can be used as in (4), 1101 . However, this method 
converges very slowly. To develop a faster algorithm, we first 
show that .^(Xp) is a non-decreasing function of /.t > 0. To this 
end, suppose /ii > /i 2 . Since R lcq — <7(X M ) is a subgradient 
of d(n) at /i, d((ii) < d(ii 2 ) + (Mi - ^2)(R Icq - sfX,,,)) 
and d(/x 2 ) < + (/*2 - A*i)(^r 0q - 5(X Ml )) follow. As a 

result, 

(Ml - M2)(^rcq - 5( X Ml)) ^ - ^(^2) 

< (Ml -^(-Rreq-ffCX^)) 

holds, and thus g(X Ml )) > ^(X^), meaning that .^(X^) is 
indeed non-decreasing with fi. Based on the above property, 
the bisection method can be used to the [i > satisfying 
c/(X M ) = R rcq as n*. 

The overall procedure to solving (P2) for X* is shown in 
Algorithm Q] where e > is small and prescribed. As can be 
shown in a similar way as in lfl6l . the finally produced X ;t 
is either equal to (if ^(Xp) = R req is satisfied), or a close 
approximation (if R req < g(X M ) < R lcq + e is satisfied) for 
X*. Moreover, the indicator variables in that X M are either 
or 1, and therefore the corresponding S is either optimum or 
approximately optimum for (PI). It can readily be shown that 
Algorithm [T] has a polynomial complexity with respect to K. 

Algorithm 1 The algorithm to solve (PI). 



Vfc,Z; 

tx = 1; compute ^(X^^ 

t ) < -Rrcq dO 

lax ; compute g(X Mmax ); 



compute Gki, 

/^min = O5 /^mt 

while g(X Mma: 

end while 
while 1 do 

^ = Mmax+Mmi„ . ^Jyg ( p 3) foj . X ^ ; 

if -Rrcq < ff(X M ) < R rcq + e then 

go to line 15; 
else if g(X M ) > i? rcq + e then 

else 



/'1 



end if 
end while 

compute the S corresponding to X M as an (at least 
approximately) optimum solution for (PI). 



is solved for its optimum solution {t kl \V fc, /}, an optimum I 
for (P4) can be constructed as follows. For every combination 
of fc and Z, the entry in {t\\ , t?, } with the smaller metric is 
assigned with t^, while the other entry with 0. 

Most interestingly, (P5) is a standard assignment problem, 
hence {^ ; |V fc, 1} can be found efficiently by the Hungarian 
algorithm, and every entry in {^ ; |V k,l} is either or 



IV. Numerical experiments 

Consider the scenario where the relay is located in the 
straight line between the source and the destination. The 
source-to-destination and source-to-relay distances are 1 km 
and d km (d G [0, 1]), respectively. The parameters are set 
as cr 2 = -50 dBm, R lcq = 100 bpos and e = 1. When K 
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and d are fixed, every channel impulse response is randomly 
generated in the same way as in |fT71 . 

To illustrate the power-reduction benefit of the improved DF 
protocol with OSP, two benchmark protocols are considered. 
The first one is the improved DF protocol with a priori 
subcarrier pairing as studied in JT]. The second one is the non- 
cooperative transmission, i.e., the direct transmission is used 
at every subcarrier. Define P sp , Pf sp and Pd as the minimum 
sum power needed for the improved DF protocol with OSP, 
the first and second benchmark protocols, respectively. Define 
JV S p and ATfgp as the optimum number of subcarrier pairs used 
for the relay-aided transmission by the improved DF protocol 
with OSP and the first benchmark protocol, respectively. P sp 
and Nm can be computed with Algorithm 1. It can readily 

be shown that Pd = 2 V, A — — , where A satisfies 

'— IR [ IJsd.fc J 

that J2k c ([^ G sd,k - = Moreover, P fsp is equal to 
the optimum objective value of (PI) imposed with the extra 
constraint = = 0, V k, I : k ^ I. An algorithm similar 
as Algorithm Q] can be designed to find Pf sp and Nf sp , which 
is omitted here due to space limitation. 
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Fig. I . The numerical results for different combinations of K and d. 
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We have computed the average P sp , Pf sp , Pd, -g£ and 
over 1000 random channel realizations for different 
combinations of K and d. The results are shown in Figure 
Q] It is shown that for any fixed combination of K and d, 
the average P sp is smaller than the average Pf sp and Pd, 
which illustrates the power-reduction benefit of the improved 
DF protocol with OSP. Moreover, the average Pf sp is smaller 
than the average Pd- This is because the first benchmark 
protocol uses opportunistic DF relaying, which better exploits 
the flexibility of transmission-mode selection for the sum- 
power reduction. 

When K is fixed, it can be seen that the average P sp and 
Pf sp reduce while the average and increase if the 
relay moves towards the middle between the source and the 
relay. This trend for the average P sp and j s explained as 
follows (the one for the average Pf sp and can be explained 
in a similar way). Obviously, the pairing of more subcarriers 
for the relay-aided transmission is more beneficial for sum- 
power reduction if V k, I, Gki is more likely to take a high 
value. Note that Gm takes a high value only if both G sr .k and 
G r d,i are much higher than G s< x,k< which can be verified by 



using the intuitive method explained in the Appendix of |8j. 
When the relay lies in the middle between the source and the 
relay, it is more likely to have G SIt k and G r d,z both be much 
greater than G s d.fe, and thus Gki is more likely to take a high 
value. This explains the observation. 

When d is fixed and K increases, it can be observed that 
the average P sp and Pf sp reduce while the average ^jp and 
— increase. Moreover, the average P sp and Pf sp are much 
smaller than Pd, and the average P sp is much smaller than 
the average Pf sp , especially when K takes a high value. This 
is because using more subcarriers leads to more flexibility 
of subcarrier pairing and transmission-mode selection for the 
sum-power reduction. 

V. Conclusion 

We have developed a sum-power minimized RA algo- 
rithm subject to a sum-rate constraint for cooperative OFDM 
transmission using the improved DF protocol with optimized 
subcarrier pairing. The power-reduction benefit of this protocol 
has been illustrated by numerical results. 



References 



[il 



L. Vandendorpe, J. Louveaux, O. Oguz et al., "Improved OFDM 
transmission with DF relaying and power allocation for a sum power 
constraint," in ISWPC, 2008, pp. 665-669. 

Y. Wang, X. Qu, T. Wu, and B. Liu, "Power allocation and subcarrier 
pairing algorithm for regenerative ofdm relay system," in IEEE Veh. 
Techn. Conf., Apr. 2007, pp. 2727-2731. 

Y. Li, W. Wang, J. Kong et al., "Power allocation and subcarrier pairing 
in OFDM-based relaying networks," in IEEE Int. Conf. Commun., May 
2008, pp. 2602-2606. 
[4] M. Hajiaghayi, M. Dong, and B. Liang, "Optimal channel assignment 
and power allocation for dual-hop multi-channel multi-user relaying," in 
Proc. IEEE INFOCOM, Apr. 2011, pp. 76 -80. 

L. Vandendorpe, R. Duran, J. Louveaux et al., "Power allocation for 
OFDM transmission with DF relaying," in IEEE Int. Conf. Commun., 
May 2008, pp. 3795-3800. 

L. Vandendorpe, J. Louveaux, O. Oguz et al., "Power allocation for im- 
proved DF relayed OFDM transmission: the individual power constraint 
case," in IEEE Int. Conf. Commun., 2009, pp. 1-6. 

, "Rate-optimized power allocation for DF-relayed OFDM trans- 
mission under sum and individual power constraints," Eurasip J. on 
Wirel. Commun. and Networking, vol. 2009. 

T. Wang and L. Vandendorpe, "WSR maximized resource allocation in 
multiple DF relays aided OFDMA downlink transmission," IEEE Trans. 
Sig. Proc., vol. 59, no. 8, pp. 3964 -3976, Aug. 2011. 

, "Sum rate maximized resource allocation in multiple DF relays 
aided OFDM transmission," IEEE J. Set Areas on Commun., vol. 29, 
no. 8, pp. 1559 -1571, Sep. 2011. 

C.-N. Hsu, P.-H. Lin, and H.-J. Su, "Joint subcarrier pairing and power 
allocation for ofdm two-hop systems," in ICC 2010, May 2010, pp. 1 
-5. 

[11] H. Boostanimehr, O. Duval, V. Bhargava, and F. Gagnon, "Selective 
subcarrier pairing and power allocation for decode-and-forward ofdm 
relay systems," in ICC 2010, May 2010, pp. 1 -5. 

G. Li, Z. Xu, C. Xiong, C. Yang, S. Zhang, Y. Chen, and S. Xu, "Energy- 
efficient wireless communications: tutorial, survey, and open issues," 
IEEE Wireless Communications, vol. 18, no. 6, pp. 28 -35, Dec. 2011. 
S. Boyd and L. Vandenberghe, Convex optimization. Cambridge 
University Press, 2004. 

[14] D. P. Bertsekas, Nonlinear programming, 2nd edition. Athena Scien- 
tific, 2003. 

H. Khun, "The hungarian method for the assignment problems," Naval 
Research Logistics Quarterly 2, pp. 83-97, 1955. 
T. Wang, F. Glineur, J. Louveaux, and L. Vandendorpe, "WSR maxi- 
mization for downlink OFDMA with subcarrier-pair based opportunistic 
DF relaying," submitted to IEEE Trans. Signal Proc, 2013. 

[17] T. Wang and L. Vandendorpe, "Iterative resource allocation for maxi- 
mizing weighted sum min-rate in downlink cellular OFDMA systems," 
IEEE Trans. Signal Process., vol. 59, no. 1, pp. 223-234, Jan. 2011. 



[21 



[3] 



[5] 



[61 



[71 



[8] 



[9] 



[10] 



[12] 



[13] 



[15] 
[16] 



